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CHAPTER 1 
INTRODUCTION 


1.1 GENERAL 

One of the most important characteristics of a nuclear 
reactor is the effect of temperature on reactivity. As the 
temperature of the moderator, fuel, or reflector are 
individually changed, the effect on the reactivity may be 
quite different. Because it is difficult to mock up a 
temperature profile in a low power reactor, it is common to 
discuss temperature effects in terms of the isothermal 
temperature coefficient of reactivity where isothermal implies 
that all parts of the reactor are at a constant and uniform 
temperature. Comparison of the measured and calculated iso- 
thermal temperature coefficient aids progress toward under- 
standing the physical processes taking place in the reactor. 

This dissertation is concerned with two aspects of 
thermal effects on reactivity. The first is an experimental 
study of the isothermal temperature coefficient of reactivity 
of the Cornell University Zero Power Reactor, a light-water 
moderated, 2.1% enriched, U0, fueled, aluminum clad reactor 
having nominal water-to-fuel ratios of 1:1, 1.5:1, 2:1, 3:1, 


and 4.1 2/ (See Appendix A for fuel and lattice parameters.) 


1 


—" §. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores,'' CURL-7, p. 220 Cornell 


University Reactor Laboratory (1964). 
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A theoretical prediction is given in Chapter 2, and the experi- 
mental techniques and results are discussed in Chapter 3. 

The second aspect of this work singled out the temperature 
dependence of the disadvantage factor as one contribution to 
the temperature dependent behavior of the reactor. Chapter 4 
discusses the theoretical basis of the use of foils to measure 
disadvantage factors. Experimental methods and results are 
discussed in Chapter 5 and compared with predictions by Meevesn =" 


Chapter 6 summarizes the results of both aspects of this 


work and suggests further experimental studies. 


1.2 ISOTHERMAL TEMPERATURE COEFFICIENT OF REACTIVITY 

The isothermal temperature coefficient a is defined as 
dop/d@ where p is excess reactivity and @ is the temperature 
of all parts of the core. J[sothermal coefficients of reactiv- 
ity have previously been measured for cores similar to the 
Cornell critical assembly. In these studies, temperature 
dependent excess reactivities were determined either by 
measuring periods for a given control rod configuration or 
by measuring the change in critical position for a calibrated 


control rod. 


Zi C. R. MacVean, "A Simplified Cell Theory Applied to the 
Calculation of Thermal Neutron Spectra in Light Water 
Lattices,'' CURL-10, p. 145 Cornell University Reactor 


Laboratory (1964). 
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The Westinghouse temperature coefficient experiment s222/ 


were made on aluminum clad, 1.3% enriched metal and oxide 
fuels of both 0.600 and 0.38/-inch diameters. Reactivity 
values were determined by measurement of periods longer than 
100 seconds. Water-to-fuel ratios from 1.5:1 to 5:1 were 
studied over a temperature range of 20-50°C. 

The Yankee cores, also studied by Weoeinehonse| had 
2.2:1 and 2.9:1 water-to-fuel ratios for 2./% enriched, 
stainless steel clad, 0.300-inch diameter U0, fuel. Reactiv- 
ity determination was by period measurement with various 
control rod configurations. The temperature range covered 
was 22-75°C. 

The Cornell University Zero Power Reactor was studied 


over a temperature range 5-70°C for water-to-fuel ratios of 


1:1, 2:1, and 3:1. After clearly establishing that the 4:1 


J. R. Brown, et al., "Reactor Properties of Water Moderated 
Slightly Enriched Uranium Lattices,'' WAPD-117, p. 117 


Westinghouse Atomic Power Division (1954). 


— J. R. Brown, et al., "Kinetic and Buckling Measurements on 
Lattices of Slightly Enriched Uranium or U0, Rods in 
Light Water,'' WAPD-176, p. 39 Westinghouse Atomic Power 


Division (1958). 


= P. W. Davison, et al., "Yankee Critical Experiments-Measure- 
ments on Lattices of Stainless Steel Clad Slightly Enriched 
Uranium Dioxide Fuel Rods in Light Water. " YAEC-94. n. 7’ 


Yankee Atomic Electric Company (1959). 
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lattice had a positive isothermal coefficient of reactivity 
with no indication of going negative, the study of this core 
was terminated at a maximum temperature of 45°C. The 
accuracy of the values of a for both the Westinghouse and 
Cornell experiments is about f 0.03 ¢/°C, and this accuracy 


is nearly independent of the value of a. 


1.3 TEMPERATURE DEPENDENT DISADVANTAGE FACTORS 

The neutron density disadvantage factor OF and the 
experimentally measured disadvantage factor O axp are defined 
by the ratios 6. = Wyn ana Da = tc! where WN and T are 
respectively the spatially averaged neutron density and foil 
count rate per unit volume; superscripts m and f indicate 
respectively the moderator and fuel regions of the reactor. 

Disadvantage factor measurements have previously been 
made on cores similar to the Cornell critical assembly. The 


ear ee ace 0) 


Brookhaven measurements———-+== were made with dysprosium 


G/ H. Kouts, "Intra-cell Flux Traverses,'' BNL-1783, 
Brookhaven National Laboratory (1953). 


7/ 


H. Kouts, "Intracell Flux Traverses and Thermal Utilization 
for 1.15% Enriched Uranium Rods in Ordinary Water," 
BNL-1982, Brookhaven National Laboratory (1954). 


8/ G. A. Price, "Thermal Utilization Measurement,'' BNL-1992, 


Brookhaven National Laboratory (1954). 
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pellets distributed through the fuel and moderator regions. 
The pellets were 1/16-inch diameter, 0.007-inch thick foils 
made from a mixture of Dy 704 with aluminum or plastic. 
Numerical integration of the detailed flux maps gave average 
count rates for each region which were converted to experi- 
mental disadvantage factors. An uncertainty of about 2% was 
assigned from data scatter; statistical uncertainty in 
counting was about 1%. 

The Westinghouse disadvantage factor Meee umenenee ae 
were made with uranium sector foils cut and placed in the 
fuel and moderator regions of the reactor. The foil was 
highly enriched y+? dispersed in aluminum to match the fuel 
ut? atom density. These foils were gamma-counted for 
fission product activity. The experimental disadvantage 
factor was obtained directly as the ratio of count rates per 


unit volume of the moderator and fuel sector foils. The 


uncertainty in the measurement was about 3%. 


9 / 


—" H. J. Kouts, et al., "Intracell Flux Traverses and Thermal 
Utilization, 1.027% Enriched Uranium Rods in Light Water," 
BNL-1796, Brookhaven National Laboratory (1954). 


aU} H. C. Honeck, "The Calculation of Thermal Utilization and 


Disadvantage Factors in Uranium/Water Lattices,'' Nuc. 


Sci. Eng. 18, 49 (1964). 
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No previously published works indicate that the tempera- 
ture dependence of the disadvantage factor has been studied 
experimentally. MacVean's recently developed method of cal- 
culating thermal lattice papameeere 2" gives the temperature 
dependence of the disadvantage factor. Since the calculated 
changes are of the order of 1077/°c an experimental method 
was required which was more accurate than those previously 
used. 

Experimental disadvantage factors were measured for 
three water-to-fuel ratios over a temperature range 5-65°C. 
An integral foil counting method was used employing a 
0.600 * 0.00l-inch diameter foil in the fuel and an appro- 
priately shaped moderator foil. Five-mil thick foils of 
dysprosium (5.24% dysprosium in aluminum), manganese (P-metal), 
and vanadium were chosen because their cross sections are 
nearly 1/v. Two counting methods were used. One was to 
dissolve the manganese foils in nitric acid and gamma-count 
aliquots of the solutions in a well-type scintillation 
counter. The other method was to beta-count the solid foils 
in a proportional flow counter. Results of the two methods 


agreed within the statistical errors which were about 0.7%. 


il/ A. Z. Kranz, "Measurements of Thermal Utilization Resonance 


Escape Probability, and Fast Fission Factor of Water 
Moderated Slightly Enriched Uranium Lattices,'' WAPD-134, 


Westinghouse Atomic Power Division (1955). 
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i2/ A. Z. Kranz and G. G. Smith, "A Second Report on Measure- 
ments of f, p, and € of Water Moderated Slightly En- 
riched Uranium Lattices,'' WAPD-151, Westinghouse Atomic 
Power Division (1956). 

13/ 


C. R. MacVean, footnote 2, pp. 143-145. 





CHAPTER 2 


THEORETICAL CALCULATION OF THE ISOTHERMAL TEMPERATURE 


COEFFICIENT OF REACTIVITY 


2.1 INTRODUCTION 

A calculation of the isothermal! temperature coefficient 
of reactivity was made using two group, two region pertur- 
bation theory. The presentation here follows closely that 
of Paecimore 

Section 2.2 presents and discusses the two group, two 
region neutron balance equations and the associated boundary 
conditions necessary for the solution. The bilinearly 
weighted. first order reactivity equations from perturbation 
theory are given for the core, reflector, and interface 
regions of the reactor. in these latter equations, tempera. 
ture derivatives of the reactor parameters are weighted by 
the fluxes, adjoint fluxes, flux gradients, and adjoint flux 
eradients. The temperature derivatives are discussed in 
Sec. 2.3. Utilization of these equations for machine cal- 


Culations 15 discussed in Sec. 2.4; 


—' J. A. Larrimore, 'Temperature Coefficients of Reactivity 
in Homogenized Thermal Nuclear Reactors," Ph.D. Thesis, 


pp. 185-211, M.I.T. (September, 1962). 
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2.2 TWO GROUP, TWO REGION EQUATIONS 
2.2-l1 Steady State Neutron Balance Equations 
The steady state neutron balance equations for the core 


region of the reactor are: 


V oe a cn — 
le ie ie ee (2.1) 
>; + (D 77a ) © = Q (Jee) 
P16%le ZC ZC Ze , : 


Analogous equations for the reflector are: 


| 2 7 
(D,.V “2 14) Ply = 0 (2) 


| De» sce — 
Zr Pir + (DoF -2 97) SQ, = 0 (2.4) 


wheres subscript 1 refers to the fast group, subscript 2 
refers to the thermal group, subscript c denotes the core, 
subscript r denotes the reflector, D is the diffusion constant 
for a given group and region, ay is the fast group removal 
cross section, Hees is the fission cross section, 29 is the 
thermal absorption cross section, ¥ is the average number 
of neutrons per thermal fission, ¢€ is the fast fission factor, 
y- is the Laplacian differential operator, and p is the 
resonance escape probability. 

The adjoint equations for the core and reflector are 
given by 


2 


=f = | | 
(D1 WV 21) Vie ? Pte. 0 (2.9) 


+ a a: Se 
VG eG lc + (DV eS * 0 (2.6) 
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2 a a ae | 
iN ae Vie ein, (2.7) 
(Dy Vi-E5.) 95, = 0 (2.8) 


where the superscript + indicates the adjoint flux functions. 
Boundary conditions at the outer extrapolated boundary 
require that the fast and slow fluxes and adjoint fluxes 
vanish. At the core-reflector interface the boundary con 
ditions require that fast and slow fluxes and adjoint fluxes 
match (i.e. Pian P19 Ve = Ole and similar equations with 


the 1 replaced by 2), and that the fast and slow currents 


D.. where i 


=F. 
? Li iis 


and adjoint currents match (i.e. D ie 


Le 
is either 1 or 2, and the same expressions are written for 


the the adjoint fluxes). 


2.2.2 Bilinearly Weighted Reactivity Relations 


Once the fast and slow fluxes and adjoint fluxes for 
the core and reflector are known, they can be used in the 
standard two group, two region, first-order perturbation 
theory reactivity equations. Perturbations of the parameters 
are equated to temperature derivatives of the parameters, 
and the reactivity effects are calculated. 


ley) 


The following expression is given by Larrimore==" for 


the core contribution to the isothermal temperature co- 


efficient, a. 


OT 


io! J. A. Larrimore, footnote 14, p. 191. 
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=] 
- : es 
oo fo V7 Ue Je / dr ue 


Cc Cc 


: Le -Vor ; 207 oer Sie 
“de V%1c “V%2c ° “de Vb2e ~ Mle “dS * 1c 


dees. 4 d(p2, =) 
o Te or Be Te de Pie 


Pee deo, 
720 “do” °2¢] (2.9) 


Mewes dr is Che volume clement," the subseript con the 
integral sign implies integration over the core of the 
reactor, and @ is temperature. 

Similarly, the reflector contribution, «a 


16/ 


by Larrimore— as 


+ ae + 
Oy = dr P16 VY Fi6 P26 dr | -Ver, 
Cc r 


os 1s given 


d> 


dD dD 
ilete + ae + Ilr 
—a~ V 9), -\/°5 + “war V 5, = P+ 67 1G 
d> d>) 
+ i Le , 
+ >, “do Pir 7 °2¢ a6 2x] : (2.10) 


16/ J. A. Larrimore, footnote 14, p. 194. 
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Finally, Larrimore gives for Os, the interface ire 
+ ee =) (7 
ies OE ie 16. 2c | -Vot, OV Die ie lie 
eC 
+ + 
Vo5y V9.0 (Dy¢-Dor) ~ 17916216721 y) 
as 
F 97910 (P™167* 1) 
+ + 
F P1rP2c Fis - "2%26 E262 9y)] 
Lie 
[dv/de] (2 
where all the functions in the second set of brackets are 
evaluated at the core-reflector interface, and dV/dO is the 
change of core volume with temperature. For an infinite 
cylinder, dV/d@ is given by 
dv/d0 = 2nr* 6.4/3 (2.12) 


where r is the core radius, and Bal is the volumetric co- 


efficient of thermal expansion of aluminum. 


i7/ J. A. Larrimore, footnote 14, p. 209. 


SOGlweL we = 265 1s 


The total isothermal temperature coefficient for the 
reactor is the sum of the core, reflector, and interface 


COnertbuULLoONS : 


a=a data. . (2515) 


2.3 TEMPERATURE DERIVATIVES OF THE NUCLEAR PARAMETERS 

In order to calculate the various expressions for the 
isothermal temperature coefficient, it is necessary to have 
relations for the temperature derivatives of the various 
nuclear parameters, viz., v,€,p, and the various D's and 
Ss 
2.3.1 Reflector Group Constants 

The temperature derivative of the macroscopic cross 


section > is easily found by differentiating 
=» = Noo , C2514) 


where N, is the number density of atoms whose microscopic 


cross section is o to give 
(1/5)d=/de = (1/N,)dN_,/de + (1/0)do/do . (Qos 


Since N, is proportional to the density d of the material, 


Lt boLtows that 


(1/N,)dN,/do = -p , (2.16) 
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where B is the volumetric coefficient of thermal expansiun 
for the material. By substitution in Eq. (2.15), the tempera- 


ture derivative of * becomes 
(1/=)d=/dO = -B+(1/0)do/do . CZ.) 
For the fast group, it was assumed that the microscopic 
cross section does not change with temperature. This 
assumption leaves only a density effect term for the fast- 


group macroscopic cross section temperature derivative. For 


the reflector we therefore have 


(1/2, ,.) dj ,/d0 = “PH0 , (2418) 


The slow microscopic cress section for water Oo, is 


proportional to 1/v; hence it is proportional to gi lle, 


Therefore, 


(1/o5,.)do5,/de = -0.5/0 (2,19) 


where 96 is in °K. For the slow group in the reflector, 


Eq. (2.17) therefore becomes. 


(1/25) di 5,/d0 = Py 772 : (2208 


The diffusion constant D can be written 
D = La) : (22 


where ey is the appropriate transport cross section. 





SeCa 25 ils. 
Logarithmic differentiation gives 
(1/D)dD/doe = +B-do,./d@ é (25225 


el 8/ gives the following relationship for the 


Larrimor 
logarithmic slope of the thermal transport cross section 


of water: 


(1/o,,.)do,./do = -n/O , (2.255 


where n is given for 20, 50, and 100°C as 0.473, 0.471, and 
0.469 respectively. Using the value for 50°C, Eq. (2.22) 


becomes 


(1/D5_.)dD5,/d0 = aggttet ise ( e2s 


Assuming again that the microscopic cross section of the 
fast group does not change with temperature, Eq. (2.22) 


for the fast group becomes 


(L/D, ,.) dD, ./de@ = PHO (2225) 


Peo. ‘Core Group Constants 


In the core region of the reactor, the fast group 


macroscopic cross section me can be written 


oe fi tie Le L Aaree 
Bee (Wve) Sets) (2.26) 


18/ J. A. Larrimore, footnote 14, pp. 83-84. 
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where ve and V" are the volumes of the fuel and moderator 
region of the cell, me is the fuel cross section, and he 
is the moderator cross section. 

For the fast group, the microscopic cross sections are 
assumed not to vary with temperature, and the expansion of 
the fuel is taken to be negligible compared to the expansion 
of the cladding. The Logarithmic temperature derivative of 


ay then becomes 
c 
= _ ai ees lie 2 ; - 
(1/2,,)42,./40 = V™view™) he /540) (264 /3) uo] > 


(2 .2f) 
where Bay is the volumetric coefficient of thermal expansion 
of aluminum, the cladding material. The last term in 
Eq. (2.27) arises because the density effect on the moderator 
is the combined expansion effects of the water and aluminum 
For the reactor cell, it is only the water cross-sectional 
area change per unit length of the reactor cell which con- 
tributes to the water density change. This area change is 
due to the expansion of aluminum structural members. The 
2B 43/3 term accounts for this area change per unit length. 


Using Eq. (2.21) where, for the fast group, 


i Goo) a oe +t st ) 


t er ee * (2.28) 


it can be shown that the logarithmic temperature derivative 


of Di. 1s 





Seey 250 ib 7 


a Tiley eee eis al eee 
Ci Diy) dD Got Dye Va) y072Pa1/3) 


(25299 
The temperature derivatives for the slow group core 
parameters were taken from calculations by MacVean using the 
computer program COUTH, "a simplified polyenergetic cell 
theory (which) is formulated to determine spatially averaged 
energy dependent thermal fluxes in the moderator, cladding, 
and fuel regions within the unit cell of a reactor lattice. 
The derived spectra are then utilized in the calculations of 
the thermal integral parameters and average cross sections 


ul9/ Calculations at tvs 


required for reactor computations. 
temperatures provide the temperature dependence of the thermal 


parameters used here 


2.3.3 Resonance Escape Probability 
The resonance escape probability temperature derivative 
can be calculated by the following scheme Weinberg and 


20/ 


Wigner——" give the following expression for p(E), the 


resonance escape probability to energy E from source energy 


sd 


C. R. MacVean, footnote 2, p. ii 


— Ae M. Weinberg and E P Wigner, Ihe Physical Theory of 


Neutron Chain Reactors, p. 301. University of Chicago 


Press, Chicago (1959) 
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E | 


5 
p(E) = exp | -(1/2) | a(E') S CE") 450] 1dE'/E! . 
E 
(20) 


where € is the average lethargy increase per collision, 2 a 
is the energy dependent absorption cross section of the cell, 
and me is the scattering cross section of the cell. Equation 


(2.31) is rewritten in the usual form 


is : 9 se! | 
P= EAP N O626) Bier (22505) 


where Rl vce is the effective resonance integral, and No we 
the homogenized fuel atom density in the cell. 
Logarithmic differentiation with respect to temperature 


results in the following expression 


= , 2 -l. ie .¢ T 
(1/p)dp/de = -N Riogst% .) i 1/N,) dN, /do (1/5. .)d ,/de + 


(L/R1 ge) IRI og ¢/ dO. 
(2.52 
The atomic density of the fuel varies with temperature. 


however, since fuel atoms are lost from the cell only 


through axial expansion of the fuel, we have 


(1/N ) dN, /do a a ae (2.233 


u 


The microscopic scattering cross section 1s a constant, 
therefore the scattering cioss section temperature deriva- 


tive is given by 
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(1/=,) dz ,/de = -B CPS) 


S 


The scatterer is water and is lost from the cell by expansion 
of water, but is gained in the cell as the horizontal cross- 
sectional area of the cell increases due to aluminum 
structural expansion. The density change in the moderator 

is therefore a combination of water and aluminum expansion, 


and 
(1/5,) d= _,/de = “Py ot (2/3)F a1 : G2 oe) 


Using Eqs. (2.33) and (2.35), Eq. (2.32) can be re- 


written as follows 


(1/p)dp/d@ = -N RI ¢e(62,)7° 


2 1 
Puno 7 3 Pal ~ 3 Pu F 
(L/RI.¢¢) dR, ¢-/ de] 
(230) 


Table 2.1 gives the values of the terms in Eq. (2.36) 


for the water-to-fuel ratios investigated. 


Z-5.4 Kast Fission Factor 

An estimate of the temperature effect on the fast 
fission factor can be made using an empirical relation 
for the fast fission factor for U0,-water lattices as 


given by areinore 


21/ J. A. Larrimore, footnote 14, p. 145. 
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(1/2,)d=_/dO = -B (2.34) 


S 


The scatterer is water and is lost from the cell by expansion 
of water, but is gained in the cell as the horizontal cross- 
sectional area of the cell increases due to aluminum 
structural expansion. The density change in the moderator 

is therefore a combination of water and aluminum expansion, 


and 
(1/5) dz 5/de = Py ot (2/3)8 41 . (2235) 


Using Eqs. (2.33) and (2.35), Eq. (2.32) can be re- 


written as follows 
(1/p)dp/de = -N RI (E> ia [s s ZB . i B + 
Pee u*eff\5~s HO 8 Ale een 


(1/RL¢¢)dR1I,--/de] 
(2.36) 


Table 2.1 gives the values of the terms in Eq. (2.36) 


for the water-to-fuel ratios investigated. 


2.3.4 Fast Fission Factor 

An estimate of the temperature effect on the fast 
fission factor can be made using an empirical relation 
for the fast fission factor for U0, -water lattices as 


given by Paeeimonet 


21/ J. A. Larrimore, footnote 14, p. 145. 
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e-1 = (0.1565) [1+0.875(dy . V"/v")+0.288 VS/vi}7> (2.37) 
2 


where € is the fast fission factor, d is the water density, 


H50 
and V> is the volume of the structural materials. Logarithmic 


differentiation, assuming vs /vt to be small, gives 
(1/e)de/de = 5.59(VE/VE) (e-1)" (6) "(2 o/3)-8,/31 (2.38) 


Using these values for the 1:l core: € = 1.03, yryyt = 120i 
@ = 20°C, & , = 2.05x10°"/°C, and 6, = 0.3x107"/°C, the 

9 u 
temperature effect on € as calculated by Eq. (2.38) is found 
to be 0.6x107°/°c. For all other cores ¢ is smaller than 
1.03 and the (ee term is dominant. The fast fission 
factor temperature dependence is therefore found to be negli- 


gible. 


255.5 Discussion of the Temperature Derivatives of the Nuclear 


Parameters 

Temperature derivatives of the nuclear parameters have two 
major contributions: one which is directly related to tempera- 
ture change (such as the variation of the microscopic cross 
sections for the slow group) and one which is indirectly re- 
lated to the temperature change by the effect which this tem- 
perature change has on the thermal expansion of the materials 
in the reactor. Over the investigated temperature range, the 
temperature derivatives of the microscopic cross sections change 


by about 10%. Over the same temperature interval, the value of 


ee os 
: 
ree. a vr / — _ a 
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Bay and B,, are constant, and dRI. ¢-¢/do changes by about 3/4. 
Therefore, there are some terms in the expression for the 
temperature coefficient of reactivity which do not change 
appreciably over the temperature interval studies. Over 
this same temperature range, the value of Pao changes from 
about 0.5x1074/°¢ to about 6.0x107*/°c. 

The above discussion suggests an expression for a of 


the form 


Qa = eee AO , (24309 


Calculated values of A and B are tabulated in Table 2.3 in 
the next subsection. The experiments described in Chapter 
3 showed the predicted linear dependence although the cal- 
culated and observed values of A and B did not agree for 


all water-to-fuel ratios. 


2,4 NUMERICAL CALCULATION OF THE ISOTHERMAL TEMPERATURE 
CORP PIGIENE OF REACTIVITY 
A computer program was developed by Ostese=” which 
solves the neutron balance equations, calculates the fluxes 


and adjoint fluxes, and calculates the isothermal tempera- 


ture coefficient of reactivity by Eqs. (2.9)-(2.11) for an 


22/ D. J. Osias, private communication (August, 1964). 
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imfinite cylinder reactor Calculations were made at 4, 20, 
40, 60, and 80°C. These calculations represent an approxi- 
mation to the first order perturbation theory in that the 
flux shapes calculated at 20°C were used to compute the 
temperature coefficient at all other temperatures. All para- 
meters of the two group two region problem change with 
temperature on the order of 10° */°c, so that changes in the 
flux shapes with temperature are expected to contribute an 
error of less than 1% to the calculated values of the iso- 
thermal temperature coefficient of reactivity. 

The input numbers required and their sources are Listed 
in Table 2.2. Results of the computer calculation are given 
in Table 2.3 for the five water-to~fuel ratios of the 


reactor. The calculated values of a are a linear function 


of PHO" 











aa 


oe 
a 


- 
7 


> 

ooo 

oer a 
= id a 


> 
ee 








wvie © 


Sey: 
a 


oO 
- 


24 


edd 2 30U}003 ‘uea APN a ) 
“Fels a yOu OOu No tot Cur 
°(€96T) ArOReATOgGeT TeuoTieN suu08 ty 


CVE ul 


TOE 


2.4 


Sec. 





/2 
/ 
/ 





"OST-S¥1 “dd “0089 “INV ~“pe puz i, “satieqgsuoD sorsAud st030vey,, “UT duel ~¢ “7. /e 
p E000 —=<ly-zeme “1010 GE 1G) OY 72 wo a 
9 gE°8 ica 0G 60° BiG De) 4-01 (ep/°ap) (°a/1) 
2 U'Tl- 9°0T- — €8°6- —L6*B= 9 L- Do /y 01 (ep /°’xp) (°&z/T) 
9 L°6- 996-6 LBB LL De /y_ot Cop/* FF xp) OM /1) 
= T8810°O0 T88T0‘O I8810°O0 TS88T0‘0 T18gto‘o uP? “br 
9 OVT'O OVT'O OFT'O OFT'O }§©O4T‘O wo *¢q 
e 670'0O 6470°0 670°0 640°O0  6%0°0 7 ue ake 
e Bo | Zo got Gor cok wd *1q 
2 S€S0‘0 7290°0 L¥l0'O €80°0 €ZT60‘0 ue re 
2 90ZZ'0 "7°0 L£Z°0 TTE'O  § OLE°0 wo %q 
q 876'0 106°0 198°0 €8°0 G9L‘°0 d 
2 Z8Z0°'0  ZS€0°O SS%0°O0 ¥260°0 96650°0 we ae 
q CLO ll “Gor ee GLO Ts eco €O°T 3 
q eu c E77 €7°7 Gig €7°7 A 
e €1v0'O O170°O0 €6€0°O S8E0°O  820°0 ue “le 
2 0°62 €°6Z 8° 6Z 0°zE O°Zy Ww a 


ese 


122 





eee lend~ Oj = dicen a 





oon 40 NOTIVINDV) YOd LAANI waLNdNod 
7 7a 








: 
/ 
a 





22 


a3ed 32x9u UO peanut jUOD 


Cle. 0 7L0°O GS0°0O Eyor i= 
ye 07 7£0°0 690° 0 Coo. 0 = 
SST*O- 7L£0° 0 LEQ; 0 OO£* 0- 
Cae} 7£0° 0 080° 0 S61’ O 
veo O- 716° O O£E8° 0 716° 0 7L£0°0 060°0 OSL°O 
cee Lo 090° 0 OL0°0 o6..1- 
607° T- 090° 0 8£0° 0 Lys iL 
£96°0- 090° 0 880° 0 it = 
Goce O- 090° 0 OOT* O cya 
eke, O- 752° 0 069s 0 woo C 090° 0 Et 0 €80°0 
[00 G- 8S0° 0 660° 0 Too 
673° 17 8S0°0 7OT° 0 OTO=c~ 
Coe l= 8S0°0 9TT* 0 9CG°T- 
ee, 0- 8S0° 0 lel © Con. 
Cyan 0= CUO.0- 6698° 0 GeO 0= 8S0° 0 O71 O SoC 0 
On > c£0° 0 Te oo 
Eas. Co Ga0. 0 Lc.’ 0 086° c- 
Gb eC = CoO 0 Lae 0 LOe.G= 
bce it CEO 0 Lo 0 oy 
nog Cu U= €76L°0 Gace O- c£0° 0 vet} 8t7° O- 


2.4 


€°¢ H1dVL 


Sec. 


lee 


Lee 


T3¢°T 








26 


2.4 


Sec. 








"OI? ‘d “TT aqouRZ003 *B82ag °-s °-s /> 
°(HO6T ‘lequeqjdas) uotjeotunumos aqyeAatid ‘setso ‘f -q /q 
) oTyQ “puejareT) ‘°09 3utyst{qng szeqqny Teotwauy ‘cz <d 
‘*pe puzy ‘sotsdyg pue AxSTuey) JO yooqpuey ‘(FeTYD-UT~103Tpy) uewspoH <q <9 /e 
SS0°0- eco 070° 0 690 = GS9°9 08 
79a 0 C200 S470°0 Calne GaGs 09 
£09°0 C200 0S0°0 OS7 0 O67 oy 
L7G GOO 950° 0 816° 0 GSO’? OC 
O7¢ 0s 8ES° T SS62°0 8ES° T E20 0 7900 cOv’ T 0 7 iy 


penutqUuods €°7 aTqe] 





CHAPTER 3 


METHODS AND RESULTS OF THE ISOTHERMAL TEMPERATURE COEFFICIENT 


EXPERIMENT 


3.1 GENERAL EXPERIMENTAL PROCEDURE 

A typical series of runs for the temperature coefficient 
experiment was made by starting at the coldest temperature and 
progressing to higher temperatures in steps of five or ten 
degrees. Smaller temperature intervals were used in the region 
where the isothermal coefficient changed sign. 

Temperature measurements were made at three times during 
a run: at startup, at approximately half power, and at maximum 
power. The average of these temperatures was assigned to the 
ris 

All control rods were at the fully withdrawn position and 
the source was removed from the core prior to collecting data 
on the stable period of the clean core. Data were collected 
over approximately three decades increase in power level. 

After the final temperature measurement, the water was 
removed from the reactor tank to the dump storage tank, heated, 
and returned to the reactor tank for the run at the next higher 


temperature. 


= Di 


..e+?- 6, ~ *9— 
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3.2 TEMPERATURE REGULATION AND MEASUREMENT 

3.2.1 Water Heating 

When the reactor is shut down, the moderator water is stored 
in the dump storage tank located just below the reactor tank. The 
reactor tank capacity is 2,000 gallons; that of the dump tank is 
2,200 gallons. A steam heater in the dump storage tank was used 
to heat the water at the rate of about one-half degree centigrade 
per minute. Since the system lacks automatic temperature control, 
it was difficult to obtain or maintain a given temperature to 


within one or two degrees. 


J¢2e0¢ Water Chilling 

The recirculation pump for the reactor demineralizer system 
circulated the water from the dump storage tank through a 
Worthington water chiller23/ (model RWW-400) and back to the 
dump storage tank. The cooling rate was approximately one-half 
degree centigrade per hour. Because of this slow cooling rate, 
a series of experimental runs was made in steps of increasing 


temperature. The lowest obtainable temperature was found tec be 


apout 5°C « 


3.2.3 Temperature Measurement 
Four twenty-gauge Leeds and Northrup (#58-1) asbestcs 


insulated, fiberglass sheathed chromel-alumel thermocouples 


placed in the core through interstitial holes in the top grid 





23/ On loan from the Department of Thermal Engineering. 
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plate were used to measure the moderator temperature. Sampiing 


of the temperature took place at radii and depth as given in 











Table 3.1. 
TABLE 3.1 
THERMOCOUPLE POSTTIONS 
Water~-to- Thermocouple Thermocouple Thermocouple Thermocoup.-« 
Fuel Ratio 
Leal eras]! les (OLY 25, ea. V 8.5-I1I TO 
depth 36 18 36 24 
2.1 radius 4./9-1] 50> Tel 95-111 24 259= 14 
depth 30 24 12 36 
Sei fradius 7.5-1V 2.5-1V (aoe 5.O-I1il 
depth 29 20 18 36 
4:1 radius 8.75-T1 c/ 2.5-1V 9.75-111 
depth 18 18 18 
a/ 


— Measured in inches from the center of the core, Roman numerals 
indicate the grid plate sector. 


b/ Measured in inches down from the top of the upper grid plate 
The horizontal midplane of the fueled region of the core is 
at 30.5 inches. 

c/ 


= This thermocouple was located in the reactor fill pipe for thi-= 
water-to-fuel ratio. 


OS A A a SL a LT A KO. 





The ice junction was located in the reactor cell next to the 
reactor, and copper multiconductor cable was used to bring the 
thermocouple voltages into the control room. A precision potent:-- 
ometer (Leeds and Northrup model 8662) was used to measure the 


voltages to t 0.001 mV which corresponds to + 0.025°C. 


- 
- 

: 
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The time required for a step temperature difference between 
the cladding and the center of the fuel rod to decay to less than 
20% was computed to be of the order of three minutes. This 
analysis was made using the standard Fourier charts for transient 
heat transfer in cylinders. Helium fills the void between the 
fuel pellet (o.d. 0.600t0.003 inches) and the cladding 
(i.d. 0.610+0.002 inches). The effect of this void was included 
in the calculations by assuming that the fuel pellet does not 
touch the cladding. This assumption results in an overestimate 
of the transient decay time inasmuch as conductive heat transfer 
from the pellet to the cladding has been neglected. The fuel 


pellet properties used in this calculation were. 


fuel density 1x6 em/em> 
thermal conductivity 0.07 watts/cm-°’C 
thermal capacity 0 063 cal/gm-°C 


the time required to take the reactor critical is several 
times the three minute interval mentioned above. About fifteen 
minutes are needed to pump the water into the reactor tank from 
the dump storage tank. The reactor tank fill pump was usually 
run for several minutes after the reactor tank was full in 
order to agitate the water in the reactor tank. Another eight 
minutes are required to raise the control rods. Therefore, 


normal startup procedure allowed sufficient time to establish 
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thermal equilibrium in the reactor tank; the maximum obsex ved 


temperature difference between thermocouples was less than 7’C. 


3 Oy) PER TOD MEASUREMENT 

Following the establishment of a good temperature profile, 
the reactor was taken on a positive period by removal of all 
control rods. With the control rods withdrawn, the reactor is 
in a clean configuration because the fuel followers in the con-~ 
trol rod elements are then at the same height as the fueled 
section of the normal fuel rods 

Reactivity errors resulting from positioning errors 
should be small because the differential control rod worth 
is a minimum at the fully withdrawn position. The reactivity 
effect of a simultaneous one-inch misalignment of all four 
control rods was found to be about 0.2¢, which is within the 
statistical uncertainty of the measurement. The specifications 
for the digital indicating system for control rod position 
call for position reproducibility to within 0.05 inch. 
Therefore, misalignment of the control rods at the fully 
withdrawn position will contribute a negligible reactivity 
error, certainly less than 0.1¢. 

Time dependent flux transients associated with control 
rod motion or source removal decay rapidly so that the actual 
reactor period is within one per cent of the stable reactor 


4 / 


period after two periods have piapced = 


Cae, ee Oe 


24? | 


Lo allow transtente 





am (2 eek Oa a ee 


J Templin, "Reactor Physics Constants,'' ANL-5800, 
2nd ed., p. 463. Argonne National Laboratory (1963). 
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to decay, the data recording system was not started until the 
power level had increased beyond that at criticality by about 
one decade, which is more than two periods in time. 

The suitable range of periods was from 25 to 100 seconds. 
Periods less than 25 seconds are too short to measure satis- 
factorily. Lack of temperature control of the water in the 
reactor tank prevented the use of periods much longer than 
100 seconds at temperatures much different from room tempera- 
ture because temperature drifts cause corresponding reactivity 
drifts. To avoid problems associated with very long or very 
short periods, the excess reactivity of the clean core was 
kept between about 7¢ and about 25¢ by adding or removing 
peripheral fuel rods. This technique was chosen rather than 
that of compensating with control rods because with the con- 
trol rods fully withdrawn, systematic errors due to control 
rod positioning and to temperature effects on control rod 
worth were avoided. 

Two BF, pulse counters (Hammer-20th Century 31EB70G) 
provided the data for the period determination. The outputs 
of these counters were recorded automatically by a Hamner 


25/ The 


Spectrometer, a combination scaler-timer-printer. 
data recording cycle was adjusted so that several data points 


were recorded during the time interval of one period. The 


25/ Daw. bebe, srOGEN@te: Me pen LO: 
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detectors themselves were positioned in the reactor instrument 


thimbles so that saturation did not occur for power levels used. 


3.4 DATA REDUCTION 
Period measurement raw data were analyzed by the computer 


code rHorit,2o/ 


a least squares fit to an exponential function. 
The computed period is then fed into a subroutine which cal- 


culates the excess reactivity by the inhour equation 


6 
ex Pere eee” >, (8, /Boge) (144.7) 77 (ae) 


i=l 
where: 
\\ is the neutron generation time 
B. is the importance weighted effective delayed neutron 
yield for the io group 
A. is the decay constant for the precursor of the gon 
group 
Perf is the effective delayed neutron fraction 


T is the stable period of the reactor 


P ox! Pere is the excess reactivity in dollars. 


Table 3.2 lists the values of the parameters used in the inhour 
equation as well as selected values of P/B age and T for the 


four water-to-fuel ratios used. 


26/ W. E. Schilling and D. H. Bryce, private communication 


(March, 1964). 
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TABLE 3.2 


REACTIVITY PARAMETERS2/ 


SS SS ae, eee fae 


W/F L:1 2:1 3:1 4-] 
B, (%) 0.0397 0.0317 0.0310 040321 
By (%) 0.2206 0.2125 0.2142 OL 2147 
B3(%) 0.2023 0.1933 0.1942 0.1946 
By (%) 0.3553 0.3955 0.3945 0.3945 
Bs (%) 0.1342 0.1229 0.1216 0.1207 
Bg (4) 0.0480 0.0440 0.0437 0.0435 
Bagg (%) 0.7943 0.8630 0.8330 0.7955 
h, (sec ”) C0224 0.0124 0.0124 070124 
ho(sec™>) 0.0305 0.0305 0.0305 0.0305 
ha(sec™™) 0.1111 0.1111 0.1111 Oanian 
M,(sec™") 0.301 O2501 0.301 0.301 
Me(sec™") 1.14 1.14 1.14 ae 
Ag(sec™*) 3.01 3.01 3.01 Sam 
N/Boge(sec) 0.004041 0.005411 0.007155 0.008976 

Reactivity ean S) for period T (in seconds) 

T(sec) 

20 0.2889 0.2799 0.2812 0.2819 
30 0.2307 0.2223 0.2234 0.2240 
4.0 0.1937 0.1859 0.1868 0.1874 
50 0.1676 0.1604 0.1612 0.1617 
60 0.1480 0.1414 0.1421 0.1425 
70 0.1327 0.1265 0.1272 0.1276 
80 0.1204 0.1146 0.1152 0.1156 
90 0.1102 0.1048 0. 1es4 0.1057 
100 0.1017 0.0967 0.0971 0.0974 
a/ So. S. Bere, footnote 1, pp. Lov, 210, 232 tecorrected tb, 


De H. Bryce). 
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The difference in reactivity for two runs was divided by the 
difference in temperature to obtain values of the isothermal 
coefficient a. The average temperature 9 and the corresponding 
value of the volumetric coefficient of expansion of ver ee 


Pouce two runs were assigned to that value of a. Bu 9 was 


2 
determined by differencing a table of temperature dependent 


27/ 


volumes of water.—— Figure B.1l is a plot of 6 vs @. Tables 


HO 
Z 
B.1 - B.4 in Appendix B give the temperature, period, excess 
reactivity, and number of fuel rods for each run and for all 
water-to-fuel ratios. Values of « and Bu 9 given in the same 
g 


tables were used in another computer code to least squares fit 


the values of a to a linear function of Bu 0° 
Zo 


3.5 EXPERIMENTAL RESULTS 

3.5.1 Temperature Dependent Excess Reactivity 

Figures 3.1 through 3.4 show the temperature dependent 
excess reactivities of the cores investigated Uncertainties 
in the values of p are within the size of the points of the 
graphs. Nominal water-to-fuel ratios are designated X:1-2 
(read X to one dash two) where X is the water relative volume 
and the -2 indicates the core configuration. All the -2 cores 


were set up in the reactor with the removable unit cell at the 


27/ C. D. Hodgman (Editor-in-Chief), Handbook of Chemistry and 


Physics, 42nd ed., p. 2143. Chemical Rubber Publishing 
Co., Cleveland, Ohio (1960). 
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center of the lattice. As a result, all -2 cores are off center 
with respect to the reactor tank. A minimum distance of more 
than two feet separated the closest fuel rods from the grid plate 
Support structure. The temperature dependent behavior of the 
reactor is not expected to be affected by the eccentricity of 


the core in the reactor tank. 


Be2<2 ILsothermal Temperature Coeliteuent of Reactivity 


The isothermal reactivity coefficient for each of the four 
water-to-fuel ratios is shown in Figure 3.6 as a function of 
temperature and in Figure 3./7 as a function of the volumetric 
coefficient of expansion of water. Figures 3.8-3.1l give these 
values of a vs PHO in greater detail for each core. The 


straight lines in Figures 3./-3.1l represent the equation 


Q = A+BB,, (322) 


90 
where A and B were determined by the method of least squares. 


Fitted values of A and B are given in Table 3.3. 


3.6 DISCUSSION OF RESULTS 

3.6.1 Dissolved Air in the Water 
Some previous temperature coefficient work was plagued 
with the formation of air bubbles on the fuel code 


28/ 
29/ 


J. R. Brown, et al., footnote 4, p. 40. 


P. W. Davison, et al., footnote 5, p. 69. 
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TABLE 3.3 
FITTED LEAST SQUARES VALUES OF A AND B FOR THE EQUATION i = 





A+B By 0 
D 
a Temperature 
W/F A(in ¢/°C) B(in 10°¢) Interval 
(AG) 
1:1 0.29t0.02 -0.435+0 .005 5 = 70 
2:1 0.5370 .02 -0.367+0.006 5 - 78 
3:1 0262 50.02 -0.268+0.005 5 - 62 
4:1 1L.2470.03 -0,138t0.0L11 7 =e 


—— ee ee ee ee eG ee 


These small air voids were thought to be responsible for some 
unexpected results in the measurement of the temperature 
coefficient. Two precautions were used by previous investi- 
gators. One was to use deaerated water to prevent the problem 
from arising. The other solution was to start at maximum 
temperature and cool the reactor. Air solubility being a 
maximum at low temperatures, by cooling the water from high 
temperatures no dissolved air would be released. 

Two precautions were taken in the present experiment to 
prevent the formation of air bubbles on the fuel rods. First, 
the fuel rods were cleaned with acetone to remove grease and 
fingerprints which might tend to collect air bubbles. Second, 


water pumping was continued after the water level reached the 
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mouth of the overflow pipe leading back to the dump storage tank 
The resulting circulation of water in the reactor tank tended to 
flush away any air bubbles (and also speed up the approach to the 
equilibrium temperature) . 

Systematic errors due to variation in dissolved air content 
of the water were estimated to be negligible. In the Cornell 
Zero Power Reactor, all water heating and cooling takes place 
external to the reactor tank. The only water temperature change 
occurring in the reactor tank is a small one (1°C or less) which 
is brought about by heating or cooling of the fuel rods and the 
reactor structure by the water as it is pumped into the reactor 
tank. 

An agitator is installed in the dump storage tank to insure 
uniform properties of the water when pumped to the reactor tank 
This mixing action aerates the water; therefore it was assumed 
that the water at each temperature was a saturated solution of 
air and water. 

"Between the temperatures of 5°C and 8°C, the density of 
water saturated with air was found to be 0.000003, em/ml less 


130/ 


than the density of air-free water. At temperatures above 


this range, the solubility of air in water is less and therefore 





30/ E. W. Washburn (Editor-in-Chief), International Grivical 


Tables of Numerical Data, Physics, Chemistry and Technology, 
lst ed., p. 26. McGraw-Hill Book Company, Inc., New York 
(1928). 


7 - 
7 
7 = 7 
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the density effect is presumably even smaller. Density changes 
due to thermal expansion are large by comparison (of the order 
Oe). 

Under the assumption that dissolved air is added or removed 
in the dump storage tank, the change in solubility of air in the 
water results in a density change which is of no consequence for 


this experiment. 


3.6.2 Temperature Dependent Excess Reactivity 


The excess reactivity has been determined for four water- 
to-fuel ratios over the temperature intervals given in Table 3 3 
Figures 3.1-3.4 give temperature-dependent excess reactivities 
for various fuel rod loadings with sufficient accuracy to plan 
further temperature dependence studies. Four thermocouples 
with a total worth of about -2¢ were in the lattices for these 
measurements. 

For all water-to-fuel ratios, the excess reactivity curves 
have an initial positive slope. The temperature at which the 
maximum excess reactivity was observed was 8.8, 15.5, and 30.2 
for the 1:l, 2:1, and 3:1 water-to-fuel ratios respectively. 
For the 4:1 core no maximum was reached; a curve of excess 
reactivity vs temperature for this core is shown in Figure 
3.11. Experimental data are shown up to about 45°C; for higher 
temperatures the plotted points are from a linear extrapolation 
of a vs 8. For an increase in isothermal temperature from 20 C 


to 100°C, the reactivity added would be 35¢. Thus, for normal 
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room temperature operating conditions, a temperature increase to 
160 © will -not.sin weselie bead Co. prompl, CFLLicality.® Since 


1/ 


the Safeguards Analysis== was based on the 4:1 core and assumed 
prompt criticality, there is no significant change in the 
analysis brought about by the determination of a positive iso- 


thermal temperature coefficient of reactivity for this core. 


3.6.3 Isothermal Temperature Coefficient of Reactivity 


The isothermal temperature coefficient of reactivity has 
been measured for water-to-fuel ratios of l:l, 2:1, 3:1, and 
4:1 over temperature intervals listed in Table 3.3. The 
accuracy of the individual measurements of © is about 0.03- 
0.06¢/°C. The accuracy determined for a from the least squares 
fitted curve is 0.02¢/°C for three lattices, and 0.03¢/°C for 
the 4:1 lattice. 

Within the accuracy of the measurements, A and B are 
different constants for each water-to-fuel ratio. Theoretical 
and experimental values of A and -B are plotted vs water-to- 
fuel ratio in Figure 3.12. The experimental curves shown 
merely connect the measured points smoothly and have no 
theoretical basis. However, they permit estimates by inter- 
polation for the 1.5:1l water-to-fuel ratio: A = 0.40 to .02 in 


¢/°c and B = -0.405+0.006 in 10°¢. 

31/ S. S. Berg, "Final Safeguards Report to the U. S. Atomic 
Energy Commission for the Cornell University Zero Power 
Reactor,'' CURL-4, p. /-3. Cornell University Reactor 
Laboratory (1962). 
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Theory predicts and experiment verifies a linear relation- 


ship between a and By 9; however, Figure 3.12 shows disagree- 


ment between the aadticed and calculated values of A and B. 
This is a well-defined experiment. All parts of the reactor 
were at the same temperature. The control rods were fully 
withdrawn, and the fuel followers provided a clean core. 
For each temperature, a stable reactor period was measured . 
There were no large perturbations introduced by experimental 
devices. The experiment does not interfere with itself, as 
do flux measurements using foils. The conclusion therefore 
must be that the theory fails to predict the proper results, 
and further refinements in the calculations are necessary. 
The only available previously published data which were 
complete enough and covered a large enough temperature range 
to permit useful plots of the isothermal temperature 
coefficient of reactivity as a function of temperature, and 
as a function of the volumetric coefficient of expansion of 
water, were for the SPERT-I reactor, a highly enriched pool 
ZY 


reactor 2% The temperature range covered was 15-97°C. 


Figure 3.13 shows both a vs 9, anda vs By 9 for the 


2 
SPERT-I reactor. The temperature coefficient of reactivity 
1s seen to be a linear function of the volumetric co- 


efficient of expansion of water for this reactor also 
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3./ CONCLUSIONS 


The excess reactivity has been measured as a function 
of temperature for water-to-fuel ratios of 1l:l, 2:l, 3:1, 
and 4:1 with sufficient accuracy to allow planning of further 
temperature dependence experiments. 

The isothermal temperature coefficient of reactivity 
has been measured for the same water-to-fuel ratios and was 
found to be a linear function of the volumetric coefficient 
of expansion of water for all lattices. Furthermore, the 
results of SPERT-I also exhibit this linearity of a vs °H,0° 
Theory disagrees with measured values of A and B, indicating 


a need for more sophisticated calculations for the isothermal 


temperature coefficient of reactivity. 





CHAPTER 4 
DISADVANTAGE FACTOR MEASUREMENT THEORY 


4.1 THEORETICAL DISADVANTAGE FACTORS 


Theoretical thermal disadvantage factors are calculated 
by taking ratios of spatially averaged thermal neutron den- 
Sities, which include neutrons of all energies below the 
thermal neutron cutoff energy Be (or below the corresponding 
velocity Vo): 

Honeek = gives the following definitions for the thermal 
neutron density, N(r); thermal neutron density disadvantage 


Pactor , OS and the thermal flux disadvantage factor, 6 


Vv 
Cc 


N(x) -| N(xr,v)dv (Aa) 
0 


p° 


(1/v") | dr N'(r) 


= mod 
6 = (4.2) 
c1/v*) [ dr N*(x) 
fuel 
5 = NNT (4.3) 
09 = OS (4.4) 
og = vo iv ee 
where 
e ve 
ead VN (v) dv/ / N(v)dv (4.6) 
0 O 
33/ 


earn. GC. Homeck, Tootnote FO. p. 50. 
-56- 
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in the fuel and moderator regions, and where: 


i Ss PEWe = pOStELon wec GOL, 
Van vi are the moderator and fuel volumes; 
nN’, wi are the spatially averaged thermal neutron densities 


in the moderator and fuel regions; 
op. oe are the spatially averaged total thermal fluxes in 


the moderator and fuel regions; and 


“e 
tn 


are the velocities averaged according to Eq. (4.6) 


in the moderator and fuel regions. 


Macveane=! used similar equations to calculate the 
thermal neutron densities, fluxes, average velocities, 
density disadvantage factors and flux disadvantage factors 
for the five water-to-fuel ratios of the Cornell University 
Zero Power Reactor. His results will be given later in 


this chapter as their use is required. 


4.2 EXPERIMENTAL DISADVANTAGE FACTORS 
The experimentally measured quantity Can which is re- 
lated to the neutron density disadvantage factor is here 
called the dysprosium, manganese, or vanadium disadvantage 
factor depending on the foils used. This quantity is defined 
by 


34/ C. R. MacVean, footnote 2, pp. 143-145. 
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count rate per unit volume in the moderator foil 


6 = 
cok count rate per unit volume in the fuel foil 
aa 
For a foil, the count rate per unit volume C is given by 
CO 
C=eas (l1/V) K(t) fe J dE > (EE) (CE, 4) 
0 
(4.8) 


where 
e is the counter efficiency; 
a is the foil self-absorption correction; 
s is the foil self-shielding correction; 
V is the foil volume (integration over r is over the 
foil volume); 
2 (E) is the energy dependent activation cross section; 
O(E,r) is the energy and position dependent flux; and 
K(t) is a time dependent function which takes into account 
the irradiation time, the time since reactor shut- 
down, and the counting time. For any given dis- 
advantage factor measurement, K(t) will be the same 
for both moderator and fuel foils and will cancel 
out of the expression for the disadvantage factor. 


In practice it need not be calculated. 
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The counter efficiency, e, will be different for the 
differently shaped foils. This effect and an experimentally 
determined correction is discussed in the next chapter. For 
purposes of this chapter, e is assumed to cancel. 

For small source-to-counter distances, the self- 
absorption factor, a, is a function of the source thickness 


35/ 


emily, Any sensitivity to geometric shape is determined 
experimentally and included in the counter efficiency 
correction factor discussed in the next chapter. For purposes 
of this chapter, then, the self-absorption factors will cancel 
in the expression for the disadvantage factor. 

The self-shielding factor, s, may vary slightly from the 
fuel to moderator foils, but this variation is neglected for 
this work and the self-shielding factors are assumed to can- 
cel in the expression for the disadvantage factor. Experi- 
mental determination of the self-shielding effect could be 
accomplished by measuring ee as a function of foil thickness 
and extrapolating to zero thickness. 

If the energy range is broken up into the thermal region 
(E varying from 0 to EA); and the fast region (E varying from 
E,. to @m), the following expression can be written for a cal- 


Cc 


culated value 6 of the experimental quantity 6 : 
calc exp 


35/ 


— WwW. J. Price, Nuclear Radiation Detection, p. 132. McGraw- 
Hill Book Company, Inc., New York, Toronto, and London 


(iy 36) 
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E CO 
€ 

eas x(ey(v™)t farl f ane.ceye" (a,x) + f a= (e)s"(E,2) 

6 = Eo 

calc BE. OO 
— Kcey(vty t fir | f dE, (E)0* (E,x) +f as, @of er | 

0 E 

C 


(4.9) 
On canceling e, a, s, K(t), and on using the following relation 


for the energy dependent, spatially averaged flux, 


3(e) = v7 faroce,x) , (4.10) 


Equation (4.9) can be written as 


E CO 
f dES_(E)¢' (E) + ; dES _(E) 9" (E) 


0 E 


_ Cc 
O ale So a eS (4.11) 


C 
f dES, (E)" (E) + f dE= _(E)® (E) 
0 E 
C 
4.2.1 1/v Thermal Cross Section Foils 
For a foil with a 1/v thermal activation cross section, 


2 CE) can be written 


x (E) = ZOVolV (aes) 


where the subscript o indicates a reference energy, and 2eVe 


is a constant. The calculated disadvantage factor then 


becomes 
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Bo CO 

f dE= vv o@(E) +f dES_(E)O' (E) 
O O a 

0 E 


Cc 
CO 


i 
Cc 
fo ocsyv sta + faz af 
0 


E, 
c 


calc, 1/v 7 


(4.13) 


Using the relations 9/v = N, 9 = C/E for E > E, (where C, 


is a constant), and 


V E 


Cc c 
f N(v)dv = f N(E)dE , (4, 1) 
0 0 


1t follows that 


@ 
—in m. -l 
ZoVON + e f dE E 2 CE) 
E 
e 


Oeale. U/y ~ OO (4.15) 


—f £ f = 
ZoVN + C, dE ES, (E) 
E 
C 


By simple manipulations it can be shown that 


Th —TT) 
1+C. RI/(v oN ) 


0 = 6 
calc,1/v n wal 
1 + C, RI/(v,o JN ) 


OW 
n 


(4.16) 


| 
OO : 
: - 
4 


a ern 
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where 
0, is defined by Eq. (Any. 
Ce, C are the epithermal flux proportionality constants; 
W is defined by the expression in the brackets; 


RI is the infinite-dilution resonance integral for 

the foil; 
v_ is the thermal reference velocity (2200 m/sec); and 
o is the activation cross section of the foil at the 


energy corresponding to Vie 


If the second term in both the denominator and numerator 


of Eq. (4.16) were small compared to l, 5 i would be 


Coe a 
equal to Oe 
The resonance integral for the foil is defined by 
OO 


RI = J Eo (E)dE . Chin 


E 
Cc 


For the experimental values of RI used in the present cal- 
culations, E, was 0.4 eV. However, the values of N, BO 
and Bo calculated by MacVean are for a cutoff energy of 
1.6 eV. Corrections for the different values of E, were 
made according to the following scheme. 
Lemans ig be given by 
L0 


Ni 6 - J dE N(E) . (4.18) 
0 
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A similar expression defines No ,: Then 
6 


1 
No.4 = Ni ¢ - J dE N(E) . (4.19) 
0.4 


levee region 0.4 < E < 1.6, let 


N(E) 


p(E)/v(E) 


C,E-'v(E)™ (42) 


where C, is selected to be that value given by EQ(E) from the 
COUTH So ecera at 0.55 eV. For v(E) in m/sec, and E in eV, 


the relationship between v and E can be written 


Oe 72 3210 "° ica) sone oes (oy 


Usime, bas. (4520) and (4.21), No Rote ent On De 


= ~6 
No.4 = Ny g-1-142x10™° ¢, (4.22) 


where N is in ome and C, 16 an hana set. 


Table 4.1 gives the values of Cea Ves and RI for the 


Enree foil materials used. 


Values of W, C., No 4, and the calculated 6 


are listed in Table 4.2 for the foils used. 


cale,1l/v 


4.2.2 Foils with a Non-1l/v Thermal Cross Section 
For a foil whose thermal cross section deviates from 1/v, 
oe can be computed from Eq. (4.11) by integrating the 


energy dependent cross sections over the spectra of the fuel 
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TABLE 4.1 
RESONANCE INTEGRALS AND 2200 M/SEC CROSS SECTIONS FOR py -OF 
Mn, and V 
ae Mn V 
RI (in barns) 4322! 14 .30/ gb/ 
ae (in barns) 26005! 13.4¢/ 4 5c! 


a/ J. L. Crandall in "Reactor Physics Constants," 2nd ed. 
(L. J. Templin, ed.) ANL-5800, p. 164, Argonne 
National Laboratory (1963). 


b/ D. J. Hughes, Pile Neutron Research, p. 139. Addison- 


Wesley Publishing Co., Inc., Cambridge, Mass., 1953. 


eh D. J. Hughes and J. A. Harvey, ''Neutron Cross Sections," 
BNL-325, pp. 7-18, Brookhaven National Laboratory 


Cig 555): 
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TABLE 4.2 


THERMAL PARAMETERS FOR CALCULATING DISADVANTAGE FACTORS AT 20°C 
FOR 1/v THERMAL CROSS SECTIONS 


W/F el Deal aa 


mo ue mo ue mo ue 
c, af Galen see) 14993934103 262174. °.7995 sO371 ee0Ze 
Ny 6 b/ Cin 107 °n/fem?) Lee 5 Lh .65 28.62 21.18 41.50) 29219 
Ny, © (in 10°n/em?) 13.67 10.77 27.68 20.27 40.54 28.27 
on 0.4 1.269 L366 1.434 
Woy 0.990 0.992 0.993 
Dyer G 1.256 1.355 Lae 
Wun 0.9353 0.962 0.967 
oun, 1/v 1.209 1.314 Lied 
Wy 0.976 0.982 0.986 
Dame e253 Leo Let 
al C, = EQ(E) @ 0.55 eV from COUTH Spectra: C. R. MacVean, 
footnote 2, pp. 136-139. 
b/ 


W. E. Schilling, private communication (August, 1964). 


c/ Calculated from Eq. (4.22). 
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and moderator. The expression for Ores can be simplified 


somewhat by defining appropriate average cross sections and 
utilizing the resonance integral information above. 


The spatially averaged total thermal neutron flux for 


both the moderator and fuel regions is given by 
E 


Cc 
o- f sma. (4.17) 
0 


Average thermal cross sections are defined by 


Bo E 


Cc 
3 -f dEo, (E)$(E) f dE@(E) (4.18) 


a 
0 0 


where o for the moderator and fuel regions are integrated 


over the respective spectra. 


It can be shown that application of the above expressions 


gives for Dale 


ot L+C RI/(o o') 
ca 
of |1+ o RI/ (ot FF) 


_ ie 
= DO al Ta) x (4.19) 
For a l/v 


where X is defined by the expression in brackets. 


thermal absorber, 
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E E 
Cc Cc 
o, és J dE seve f dE (E) 
0 0 
7 Ovo /V (4.20) 


and Eq. (4.19) reduces to Eq. (4.16) as it should. 

Both manganese and vanadium thermal activation cross 
sections are very nearly 1/v. The thermal activation cross 
section for dysprosium, however, does vary from 1/v. The 
dysprosium cross sections were taken from Sher, et aijel! 
These dysprosium cross sections and 1/v cross sections for 
manganese and vanadium were averaged over the COUTH spectra 
for water-to-fuel ratios of 1l:l, 2:1, and 3:1. The values 
of these averaged cross sections, values of X, 9.4 (total 


thermal flux for a cutoff energy 0.4 eV), and the various 


disadvantage factors are given in Table 4.3. 


4.3 TEMPERATURE DEPENDENCE OF THE DISADVANTAGE FACTOR 
MacVean calculated the thermal neutron density and 
flux disadvantage factors according to the equations of 


Section 4.1. By calculating these thermal properties of 


37/ 164 


R. Sher, et al., "Low Energy Cross Sections of Dy 


Nuc. Sel. Ene. Il, 367-376 (1901). 
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TABLE 4.3 
THERMAL PARAMETERS FOR CALCULATING DISADVANTAGE FACTORS AT 20°C 
FOR NON-1/v THERMAL CROSS SECTIONS 
W/F EB Doli 3:1 


moO ue moO ue moO ue 
BN, ¢ 2/ (n/em3) = 14.58 11.65 28.62 21.18 41.50 29.19 
V1 6 2! (m/sec) 3667 3971 3134 3446 2945 3265 
31.6 (n/em?-sec) 5.3465 4.6262 8.9695 7.2986 12.2218 9.5305 


c, a (avem>—see) 0.7993 0.7708 0.8274 0.7955 0.8371 0.8028 


i a! (paens) 1735 1610 1937 1765 2019 teeo 
S, > 2! (barns) 9.300 8.733 10.22 9.452 10.60 9.711 
aie 2’ (barns) 3.507 3.295. 3<856 34566 3.999 mee: 


ay 
Fo .4 (n/em*-sec) 4.2764 3.5943 7.8618 6.2336 11.1011 8.4558 


On On4 Eel90 1,261 1.3 
Xpy 0.989 02991 0.993 
ODy 1.269 Pore 1.446 
Xun V3 0.961 0.967 
On 12208 Lo10 L206 
Xy 0.980 O953 0.98/ 
Oy 1.241 1.340 1.414 


a/ 


W. E. Schilling, private communication (August, 1964). 


bi C. R. MacVean, footnote 2, p. 141. 
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the reactor at two different temperatures (20° and 40°C), the 


temperature dependent behavior of these properties was pre- 


dicted. Table 4.4 gives the calculated values of the neutron 


density disadvantage factor at 20°C, the logarithmic slope 


(1/5) (d65_/de), and the slope db _/de. 


TABLE 4.4 


TEMPERATURE DEPENDENCE OF OF 


W/F 5 a/ (1/8_,) (d6_/do) 2! 
in 1074 /°¢ 

1:1 Lae ~2.76 

2:1 1.351 -4 .46 

ae 1422 -5.56 


a/ C. R. MacVean, footnote 2, pp. 141-143. 


do_/do 
a 10° TPG 
eG 


-6.02 


-/.91 
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EXPERIMENTAL METHODS AND RESULTS FOR THERMAL DISADVANTAGE 
FACTORS 
5.1 REACTOR DESCRIPTION 

5.1.1 The Unit Cell 

The lattices of the Cornell University Zero Power 
Reactor consist of an upper fuel rod support grid plate, 
a lower fuel rod positioning grid plate, structural support 
members, the reactor tank, the fuel rods, and four control 
rod units each containing a cluster of three control ele- 
ments in an equilateral triangle of the same pitch as the 
rest of the core. These control rod clusters may be 
positioned in any location of the reactor because the 
control elements have the same outside diameter as a normal 
fuel rod. The upper portion of each control element con- 
tains the boron carbide poison. The lower portion is a fuel 
follower with the same dimensions as the normal fuel rods. 

The upper grid plate has an insert which allows the 
removal or insertion of four fuel rods as a unit in what 
is called the unit cell. These four fuel rods form a 
parallelogram of two equilateral triangles of the lattice 
pitch. For all water-to-fuel ratios, the unit cell is 
slightly off-center with respect to the upper grid plate 


and the reactor tank. All cores studied were assembled 
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with the unit cell at the center of the fuel loading to reduce 
the radial flux gradient in the unit cell. With a triangular 
(or hexagonal) lattice, one fuel rod must be the exact center 
of the core. Therefore, of the four fuel rods of the unit 

cell, one was the central element of the reactor and the other 


three were in the first hexagon out. 


5.1.2 Activation Rods 

The original design of the reactor facility included de- 
cappable fuel rods for the purpose of introducing foils into 
the fuel region of the reactor. These rods were slightly 
longer than the normal fuel rods but contained a normal 48- 
inch fuel loading. 

Because the full-length decappable fuel rod is too long 
to fit into the hood used for handling radioactive materials, 
and because it is difficult to insert and remove foils with 
this length of fuel rod, a special "activation" rod was de- 
Signed. The eighteen-inch fuel section of the activation rod 
is supported by screw-in aluminum rods which make this special 
rod the same over-all length as a standard fuel element. 

(See Figure 5.1) When detached, the fuel section is short 
enough to be handled easily in the hood. It is easy to insert 
and remove foils with this rod because only a few fuel pellets 
need be handled. Because of the small size, ease in handling, 
and small number of fuel pellets handled, experiments can be 
accomplished with less radiation exposure to personnel than 


was previously possible with the longer rods. 





=o = 


| ___HORIZONTAL MIDPLANE OF FUEL IN 
+ F"~ THE. REACTOR 
A 


BOTTOM FUEL PELLET IN THE 
ACTIVATION ROD 


A IS THE DISTANCE DOWN FRCM THE 
REACTOR HORIZONTAL FUEL MIDPLANE 
TO THE BOTTOM FACE OF THE BOTTOM 
FUEL PELLET IN THE ACTIVATION RODS. 


Figure 5.1: The Activation Rod. 
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The activation rods are used for foil and wire activation 
within the fuel region of the core. They can be placed in any 
normal fuel rod position, but for the disadvantage factor work 
they were used only in the unit cell. The four fuel rod 
positions of the unit cell were filled with three activation 
rods and one normal fuel rod. The standard loading was as 


shown in Figure 5.2. 


5.2 FOIL PREPARATION 
59.2.1 Foil Types and Manufacture 

The foils used for the disadvantage factor measurement 
were five mil-thick dysprosium, manganese, and vanadium. 

The dysprosium was obtained from General Astrometals Corpora- 
tion and consisted of 5.24% dysprosium in aluminum. The 
manganese was in the form of an alloy called P-metal which 

is about /2% manganese, 18% copper, and 10% nickel. The 
P-metal was obtained from Metals and Controls, Inc., and has 
a thickness tolerance of +0.00025 inch. The vanadium was 
obtained from H. Cross Co. 

The fuel region foils were punched from stock with a 
0.600+0.001l-inch diameter die. Moderator region foils were 
punched from stock using a machined die of appropriate 
Shape. Figure 5.3 shows the dimensions and shape of the 


moderator foils. 
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A ACTIVATION ROD IN CENTRAL FUEL POSITION 
OF THE UNIT CELL. 


B,C ACTIVATION RODS IN FIRST HEXAGON POSITIONS 
OF THE UNIT CELL. 


D NORMAL FUEL ELEMENT IN A FIRST HEXAGON 
POSITION OF THE UNIT CELL. 


E REMOVABLE UNIT CELL FUEL ROD SUPPORT 
PLATE. 


F NORMAL FUEL RODS IN NORMAL LATTICE POSITION. 


Figure 5.2: Top View of the Unit Cell. 
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Figure 5.3: Moderator Foil Dimensions. 
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5.2.2 Foil Placement in the Reactor 

To prevent fission product contamination of the fuel 
foils, a special fuel pellet with an aluminum foil glued 
to one end with G.E. Glyptal was placed on each side of the 
fuel foils. See Figure 5.4. Ten fuel pellets were removed 
from the activation rod and an aluminum lined fuel pellet 
was inserted. The manganese fuel foil was placed next to 
the fuel pellet liner. Pellet and foil were pushed to full 
insertion with a lucite rod which was marked for trans- 
ferring the foil position to the outside of the activation 
rod. Another lined pellet was inserted, followed by a 
normal fuel pellet. This loading sequence was continued 
for dysprosium and vanadium foils until each activation 
rod had a manganese foil at the center of the rod, a 
dysprosium foil three inches below the center of the rod, 
and a vanadium foil another three inches lower. Figure 5.1 
gives the distance between the horizontal midplane of the 
fuel in the reactor and the bottom edge of the bottom fuel 
pellet in the activation rod for the three water-to-fuel 
ratios studied. 

Two moderator foils were secured to a one-eighth-inch 
thick lucite moderator support plate with double-backed 
Scotch tape. The unit cell was assembled with three acti- 
vation rods and a normal fuel rod. Using the marked lucite 
fuel foil insertion rod, the lucite moderator foil supports 


were positioned at the proper levels on the outside of the 
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Figure 5.4: Fuel Foil Placement in the Activation Rod. 
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activation rods. The assembly was taped together with Band 
Aid Air Vent adhesive tape placed about three inches above 
the manganese foils and about three inches lower than the 
vanadium foils. (This tape was found to be easy to use, to 
keep its adhesive properties under water, and to come cleanly 
off the fuel rods leaving no sticky residue.) To insure 
positive positioning of the lucite moderator support plates, 
and to prevent loss of foils, waxed string was looped 

around the fuel rods immediately above and below the support 
unit. Figure 5.5 shows the assembly ready for insertion 
into the reactor. 

The assembled unit cell was then loaded into the reactor 
and irradiated for a period of about ten to twenty minutes 
at flux levels corresponding to about one-half watt in power. 
Irradiation times and power levels were adjusted to give 
good count rates with minimum personnel exposure during 
handling of the fuel. 

It was necessary to remove the unit cell assembly as 
rapidly as possible to get any information at all from the 
short-lived vanadium foils (half life of 3.8 minutes). 

After a wait of about five minutes from shutdown, the acti- 
vation rods were removed to the radioactive handling hood 
for disassembly. The vanadium foils were removed, washed 
in acetone, and beta-counted on an automatic sample changer. 


About one hour after shutdown, the manganese and dysprosium 





Oe 
OP EDGE OF TOP FUEL PELLET 


ie _ ow 


ADHESIVE TAPE 






9" | 
Les Sua 
Ccny7 FOILS ~~~ 
FIFI IT sd ois 
3 a8 
MODERATOR FOILS 
Dy ge \rous |, 
TTS (aay a 
x 
FOILS | 
ee fa Tani 


LUCITE 
MODERATOR FOIL 
SUPPORT PLATE 


a 
ADHESIVE TAPE _ 
a i 
| L_ 


—_.— BOTTOM EDGE OF BOTTOM FUEL PELLET- 





Figure 5.5: The Unit Cell Foil Assembly. 
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foils were removed from the fuel, washed in acetone, and 


counted on the automatic sample changer. 


3-3 COUNTING TECHNIQUES -- SOLID FOILS 


5.3.1 The Automatic Counter System 

The automatic counting system described by Rosh dae 
and Bere! was used for the solid foil counting. There 
are 24 sample positions on a round table which pneumatically 
advances on electronic signal. The samples are held on 
pedestals which are pushed up next to a thin-window pro- 
portional flow counter for counting. The timing interval, 
table rotation, two scalers, and a printer are controlled 
electronically. 

Uniform counting geometry for the fuel foils was pro- 
vided by pedestals which were machined to accept 0.600- 


inch foils. The moderator foils were centered by eye on 


Similar pedestals. 


5.3.2 Relative Fuel and Moderator Foil Counter Efficiencies 


Counter efficiency is slightly sensitive to distance 


from the center of the pedeeeaia 2! Therefore an experi- 

38/ J. C. Rushing, "Report of Work Done on an Automatic 
Sample Changer for Use in the Cornell University 
Nuclear Reactor Laboratory,'' Department of Engineering 
Physics Senior Project (1962). 

oy 


—' S. S. Berg, footnote l, pp. 101-10/. 
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mental geometric correction factor was determined for this 
effect as follows. Fuel and moderator foils were irradiated 
on a flux wheel in the uniform thermal flux of the TRIGA 


wanientae! 


and both shapes were counted on the automatic 
counting system. After irradiation in a uniform flux, the 
two shapes of foils should have the same activation per unit 
volume. The foils were counted using the Geneva motion 
described in Section 5.3.4. The ratio of fuel foil count 
rate per unit volume to the moderator foil count rate per 
unit volume is the desired geometrical correction factor. 


These correction factors are given in Table 5.1 for the three 


water-to-fuel ratios studied. 


TABLE 5.1 
EXPERIMENTALLY DETERMINED COUNTER EFFICIENCY CORRECTION 
FACTORS 
W/F 1:1 2:1 3:1 
Peo! 0.990 1.007 1.014 


a/ The statistical error is about 0.004 for all factors. 


40/ E. B. Fehr, internal memo RM-16, Cornell University 


Nuclear Reactor Laboratory (1964). 


41/ E. B. Fehr and E. P. Whitted, internal memo RM-L5, 


Cornell University Reactor Laboratory (1964). 
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5.3.3 Data Reduction for Vanadium Foils 

The short half-life of the vanadium foils made it necessary 
to analyze the data differently than the data from the other two 
foils. The vanadium foils were counted in a given order on both 
proportional counters of the automatic counting system. The 
count rate at the start of counting was about 1,000 counts per 
second. The foils were counted down to about 20 counts per 
second. SAT 1, 42/ a local computer code, was used to correct 
each datum point back to the activity on removal from the 
reactor. These count rates on removal were then converted to 
count rate per unit volume. The moderator foil count rates 
were averaged, as were the fuel foil count rates. The ratio 
of these averaged count rates per unit volume for the moderator 
and fuel foils was corrected by the appropriate geometrical 


correction factor of Table 5.1; the result was the desired 


vanadium disadvantage factor Op. 


9.3.4 Data Reduction for Manganese and Dysprosium Foils 
Manganese and dysprosium have half-lives of 2.59 hours 
and 2.42 hours respectively. This is sufficiently long to 
allow use of the Geneva motion counting technique. In this 
method, the foils are counted in a given order once through. 
The last foil is counted twice. The remaining foils are then 


counted in the reversed order. For a total counting time 


42/ W. E. Schilling, internal memo, Cornell University Reactor 


Laboratory (1962). 
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sufficiently short compared with the half-life of the foils, 
a simple arithmetic time average of two counts is very close 
to the actual counts at the midpoint in time. For example, 
for a total counting time of 20 minutes, the arithmetic 
average of the end point counts is within 0.3% of the true 
counts at the midpoint time. The total counting time was 
usually about 10 minutes; therefore the error is even less. 
Statistical accuracy of counting was generally of the order 
0.7% and thus the error introduced by using the Geneva 


counting method is much less than the statistical error. 


5-4 COUNTING TECHNIQUE -- DISSOLVED FOILS 

The manganese foils were counted with a "soluble" 
method in which the foils were dissolved in acid and 
aliquots of the solution were gamma-counted in a well- 
type scintillation counter. Counting of aliquots of 
solutions of foils provides two advantages over solid 
foil counting. First, the moderator and fuel foils are 
counted under identical counting geometry. Thus, there 
is no need to determine a counter efficiency correction 
factor. Second, the process of dissolving the foil 
averages the count rate of the foil over its volume. The 
count rate of the aliquot is then directly proportional 
to the total count rate of the foil itself, the constant 
of proportionality being the fraction of the total solution 


which is being counted. For both moderator and fuel foils, 


>> = 


oa 
ae a no a 
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the constant of proportionality was the same. The above total 
count rate, when divided by the volume of the foil, is exactly 
the spatially averaged count rate per unit volume needed to 


determine the manganese disadvantage factor. 


5.4.1 Foil Preparation 

After the solid foils had been counted as described in 
Section 5.3.4, they were removed to the fume hood in the 
isotope handling room. Each foil was cut with metal shears 
so that it would fit down the neck of a ten-milliter 
volumetric flask into which about five milliliters of 
3N HNO 3 had previously been placed. Once the foil was dis- 
solved, the solution was diluted to 10 milliliters with dis- 
tilled water. Each solution was then mixed well. 

Three milliliters of the solution were pipetted off 
and placed into a plastic vial, and the vial was capped. 
Each foil had its own vial; each vial had three milliliters 
of solution. Two pipettes were used, one for the fuel 
foils and one for the moderator foils to prevent mixing in 


the pipette of solutions with different activities. 


5.4.2 Counting of the Solutions 

The vials fit snugly into the well of the scintillation 
counter, the top surface of the solution being below the top 
of the crystal. The geometrical arrangement in the counter 


was identical for each solution. 
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Three counting methods were initially investigated. 
First, a multichannel analyzer was used to get the gamma 
spectrum of the solutions. The area under the peak was 
determined by numerical integration. These total counts 
were divided by the foil weight to get counts per unit 
weight, and the ratio of the average counts per unit 
weight for moderator and fuel foils gave the experimental 
disadvantage factor. 

The second method was to count the solution on a 
Single channel analyzer in the differential mode, adjusting 
the base line and window width to straddle the gamma peak. 
Again, the ratio of counts per unit weight for moderator 
and fuel foils was the desired disadvantage factor. 

The third method used an integral discriminator with 
the base line in the valley above the Compton edge. In 
this manner, all counts above the base line were included 
in the count rate per unit weight. 

All three methods gave the same result within experi- 
mental error. The third method was chosen because of 
Simplicity of analysis, greater availability of the 
associated electronic equipment, and higher count rate 
which resulted in greater statistical accuracy. 

The aliquots were counted using the Geneva motion pre- 
viously described. Typical counting times were one minute 
per vial; typical total counts were of the order 50,000. 


The background counts were subtracted from the total counts 
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prior to averaging the counts per unit weight for the fuel 
and moderator foils. A typical data reduction is given in 


Pape’ 5.2. 


9.5 EXPERIMENTAL RESULTS 
Disadvantage factors were measured as previously de- 
scribed at various temperatures for three water-to-fuel 
ratios. The experimental results are given in Table 5.3 
All 


ag 0 1 (manganese soluble), and 5 


Dy’ Mn? °Mn, so Vv" 
values except the dissolved foil experiment have been 
corrected for geometry by the correction factors of Table 
9.L. Soluble manganese experiments had identical counting 
geometry and therefore needed no geometric correction 
factor. 

Values for the disadvantage factors from Table 5.3 
were least-squares fitted to a linear function of tempera- 
ture. Table 5.4 gives fitted experimental values and 
theoretical values of the disadvantage factors at 20°C. 
Both experimental and theoretical values of the slope of 
the 5 vs 9 curves are also given. 
and 0 


Temperature dependent values of 6 are 


Mn 


plotted in Figure 5.6 for the three water-to-fuel ratios 


Mn,sol 


investigated. Figure 5./ gives 6 vs temperature for the 


Dy 


3:1-2 core. Figure 5.8 gives 6,, vs temperature for the 


V 


same COre. 
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DATA ANALYSIS FOR A TYPICAL MANGANESE DISADVANTAGE FACTOR 


USING DISSOLVED FOIL TECHNIQUE MEASUREMENT 


Single Channel Analyzer Settings, Integral Mode 


Fine Gain Oz Base Line Li oy 
Coarse Gain 4 Window Width Oey 
Clipping Time 0.5 (microseconds) Discriminator 8.5 V 


/ Counts/ Average Average 

Foil= Foil Counts Unde Counts/ Counts/ 

Number Weight per Weight/ Weight, Weight, 
(gm) Minute Minute Moderator Fuel 


°Mn, Sol 


Foils Foils 

119 Ogl7ao 21402 237 719 310 868 236. 376 Leooas 

307 O7533 3989/ 308 366 0.007 

124 Orly43g Zia 239 273 

S17 027534 40077 3153) 370 

131 OVlI749 21375 238 136 

3d 20629 

oie 39083 

124 20774 

30 / 38564 

L119 Z2052/ 

a/ Foils numbered in the 100's are fuel foils; those in the 

300s are moderator folls. These £o1is are for the 
3:1 core. 


b/ 


— Counts per unit weight per minute are the sum of counts 
for each foil, corrected for background of 1/6 counts 


per minute, and divided by the foil weight. 
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Figure 5.7: Dysprosium Disadvantage Factors vs Temperature; 3:1-2. 
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Figure 5.8: Vanadium Disadvantage Factors vs Temperature; 3:l-2. 
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5-6 DISCUSSION OF RESULTS 
5.6.1 Foil Construction and Positioning Errors 

Moderator foils were punched from stock using a machined 
die of appropriate shape. The radius of the circular cuts 
taken in the corners of the triangles was that of a standard 
drill size (0.68/75 inch) rather than the exact radius of the 
outside of the fuel rods (0.666 inch). This gave all the 
moderator foils a slightly smaller volume than the actual 
moderator volume of interest. 

An estimate of the error associated with inexactly cut 
foils can be made as follows. First it is assumed that the 
foil thickness is uniform. A second assumption is that the 
only error in foil construction will be in making the modera- 
tor foil. This is justified somewhat in that the fuel foil 
is made from a die with less than 0,5% tolerance. The 
following equation is written for the correct experimental 


disadvantage factor for a foil with 1/v absorption cross 


section: 
—in ™m the ™m 
7 iS — MN go41% e011 + SneYne)/Y ie 
correct et cit f it ; 


(Negi foi /V 


where the subscript "foil'' indicates a property of the foil, 
and the subscript nf indicates a property of the missing 
portion of the foil, or the "non-foil.' All other symbols 


have been previously defined. The second assumption above 





Sec. 5.6 94 


sets o = Vou Equation (5 1) can then be rewritten as 


2 2 m m ym ym om 
“correct meas Poa /Y ) + (Ne at eon D Ore 


where oN. is the previously reported disadvantage factor 
as given in Table 5 3 Assuming that the missing section 
of the foil is from the circular region near the fuel 
cladding, the neutron density cladding disadvantage factors 


43/ 


of MacVean——’ can be used to approximate N./N-_.,. This 
He 101 | 
will result in an overestimate of the correction factor. 


Let 


R = v'/Ve Vfoil (3%) 


Then, Eq. (5.2) becomes 


x g ~|- 4 
“eorrect — meas Es 1L4(R-1)/9 01 (5 4) 


where CO is the cladding disadvantage factor. 

Typical values of We al teat were taken from foils 
used in the experiments. Table 5 5 gives the values of 
the corrections for the three foils and three water-to- 
fuel ratios used. [It can be seen that the correction is 
small, and generally within the accuracy of the experiment. 
Because it is a small correction, and because it is not an 
accurately determined correction, it has not been applied 
to the experimentally measured disadvantage factors 


43/ 





C. R MacVean, footnote 2, p. 141 
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Errors resulting from improper vertical positioning of 
the moderator foils in the experimental assembly are 
expected to be negligible. The flux change in the axial 
direction at the foil position in the unit cell is of the 
order of one or two per cent per inch. Positioning accuracy 
is estimated to be about 0,1 inch. The resulting error is 
then seen to be about 0.2% in count rate. 

A small error may be introduced by the presence of 
the lucite support plate for the moderator foil. This 
plate is one-eighth-inch thick and of approximately the 
same shape as the moderator foils. It is suggested that 
further experimental work be done with different thicknesses 
of lucite to investigate the magnitude of this effect. 

The error associated with the finite thickness of the 
foils has not been investigated experimentally. It was 
assumed in Chapter 4 that the self-shielding factors for 
the fuel and moderator foils were equal. In conjunction 
with the investigation of lucite thickness effect, it is 
suggested that the disadvantage factor be measured as a 
function of detector foil thickness. 

The process of dissolving the foils for the "soluble" 
disadvantage factor measurement is not expected to con- 
tribute any error to the measurement. The chemical handling 
technique employed standard volumetric flasks and volumetric 
pipettes. The pipettes are estimated to be accurate to 


about 0.25% and the flasks to about 0.1%. The error should 
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be uniform for each aliquot, and therefore cancel from the 
ratio which is the disadvantage factor. 

Weighing of the foils was accurate to 0.1 mg which 
represents a maximum error of less than 0.3% for moderator 
foils of the 1:1 core. For all other foils the error is 


tess than 0.27. 


5.6.2 Foil Induced Flux Perturbations 

There are two perturbations associated with inserting 
foils into the reactor for measurement. One effect is the 
depression of the flux by the foil itself as in the case 
of an absorber being placed into the low absorption region 
of the moderator. This would tend to reduce the measured 
value of the flux in the moderator. This effect could be 
determined by measuring the foil activity per unit weight 
as a function of foil thickness and extrapolating to zero 
foil thickness. 

The second effect is present in the fuel foil. Here, 
a highly absorbing region is replaced by a lightly absorbing 
foil. Neutrons may stream into the foil and the average 
flux in the foil would then be higher than that in a 
corresponding fuel region. This effect could also be deter- 
mined by extrapolation to zero thickness foils. 

Both of these effects tend to make the measured dis- 
advantage factor lower than the ratio of the neutron den- 


Sities appropriate to the two regions under investigation. 
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An estimate of the first correction, called the 
outer flux perturbation, was made using the procedure 


44 / 


outlined in ANL 5800. A correction factor, defined 
as the ratio of the activation in the perturbed flux 
to activation in the unperturbed flux was calculated 
for the three foil materials and for water-to-fuel 
ratios of 1:1 and 3:1 using the approximations given 
to the transport theory calculations of Ritchie and 
Eldridge. This correction is for a coin shaped foil, 
but the moderator foils used here were of the shape 
shown in Figure 5.3. An equivalent radius for these 
foils was taken to be 0.5 cm and 1.0 cm for the l1:l 
and 3:1 cores respectively. The correction was then 
estimated to be: 0.3-0.8% for dysprosium, 0.6-1.6% 
for manganese, and 0.1-0.4% for vanadium where the 
lower values are for the 1:l core. The correction 
increases the measured disadvantage factors by this 
amount. 

Because of the smallness of the corrections, and 
the uncertainty in their accuracy, they have not been 
applied to the measured disadvantage factors. [It is 
recommended, rather, that an attempt be made to measure 


them experimentally. 


B4/ L. J. Templin, footnote 24, p. 6/0. 
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5.6.3 Flux Perturbations Due to the Activation Rods 

Axial flux scans taken on the 1:1 core showed a large 
flux peaking near the ends of the fueled section of the 
activation rods. This peak is a result of replacing the 
fuel region of a normal fuel rod with the aluminum support 
pieces of the activation rod. Peaking was most pronounced 
in the 1:l core because of undermoderation for this water- 
to-fuel ratio. The peaking may have been sufficient to 
affect the vanadium and possibly the dysprosium disadvantage 
factors in the 1:l core. The manganese disadvantage factors 
are not expected to be influenced by this peaking because 
these foils were in the center of the fueled section of the 
activation rod. 

Peaking was less pronounced in the 2:1 and 3:1 cores. 
In addition, the fueled section of the activation rods was 
lowered in successive runs as indicated in Figure 5.1. 
The first change was made so that the vanadium foils were 
in the center of the core. This was changed again to the 
"symmetric" activation rod wherein the center of the fueled 
section of the activation rod is at the center of the fuel 
region of the core. 

The effect of the peaking may result in some inaccuracy 
in the measured value of the disadvantage factors for 
vanadium and possibly dysprosium, but the temperature de- 


pendence of the disadvantage factors should not be affected. 
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The only way to remove this effect is to resume use of the 
full length decappable fuel rods. A minimization of the 
effect can be achieved by using only one set of foils and 
placing them at the center of the fueled region of the 


activation rod. 


5.6.4 Comparison of Experiment and Theory 


Table 5.5 compares theory and experiment and shows 
that (1) dysprosium disadvantage factors lie between 4 
and 8% lower than theory, (2) manganese disadvantage 
factors lie between 2.4 and 4.5% lower than theory, and 
(3) vanadium disadvantage factors lie about 4% lower than 
theory. However, all measured values of the temperature 
dependence of the disadvantage factors agree with the 


theoretical values within the limits of experimental error. 


+3/ CONCLUSIONS 

The disadvantage factor was measured for three 
water-to-fuel ratios using foils of three different 
materials. Values of the measured disadvantage factors 
are consistently a few per cent lower than the calculated 
values. Further refinements in the experimental technique, 
specifically utilization of different thickness foils and 
moderator support plates, could reasonably be expected to 
bring the experimental values of the disadvantage factor 


into closer agreement with the calculated values. 
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The temperature dependence of the disadvantage factor 
was measured for three water-to-fuel ratios. Measured 
values of the temperature dependence agree within experi- 


mental errors with the theoretical values as calculated by 


45/ 


MacVean .— 


o3/ C. R. MacVean, footnote 2, p. 143 





CHAPTER 6 


SUMMARY AND CONCLUSIONS 

6.1 ISOTHERMAL TEMPERATURE COEFFICIENT OF REACTIVITY 
6.1.1 Summary 

The isothermal temperature coefficient a is defined 
as do/dO where »~ is reactivity and @ is the temperature 
of all parts of the core. The isothermal temperature co- 
efficient was measured by differencing, at a series of 
closely spaced values of 6, values of p computed from 
observed clean-core periods using a local code RHOFIT .22/ 
Temperatures were determined by averaging readings from 
four thermocouples distributed through the moderator; 
auxiliary tests showed the temperature could be made 
Ui born whem 20.) Cs 

A preliminary calculation of @ was carried out 
using two group, two region, first order perturbation 


6 / 


theory as outlined by Paerinecei= A computer program 
calculated the fluxes and adjoint fluxes and used them 
in the standard perturbation theory reactivity equations 
to find the isothermal temperature coefficient. 
Dependence of a on temperature can be ascribed to two 


contributions: a rapidly varying one proportional to 


the volumetric coefficient of thermal expansion of water, 


45/ 


W. E. Schilling and D. H. Bryce, private communication 
(March, 1964). 
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and a second one, nearly constant over the investigated 
temperature range, which is due to temperature dependence 

of microscopic cross sections and other small temperature 
dependent effects. These considerations lead to an equation 


for a of the form 


gS said: : (62 1y) 


where Py 0 is the volumetric coefficient of thermal expan- 
Sion for water. 

Experimental temperature coefficients were determined 
for nominal water-to-fuel ratios of 1:1, 2:1, 3:1, and 
4:1. Statistical errors in values of a are generally less 
than 10%. For the cores studied, and over the temperature 
eances. used (lil, 5-70°C: 2:1; 5=/8°C;: 3:1, 5-62°C; 

4:1, 7-45°C), the isothermal temperature coefficient is 
well represented by Eq. (6.1) where A and B are different 
constants for each core. A comparison of theoretical 
values of A and B with experimental values found from a 
least-squares fit to Eq. (6.1) shows a need for more 
sophisticated calculations to obtain better agreement with 


experiment. 


6.1.2 Conclusions 
Temperature dependent excess reactivities for various 


fuel loadings were determined for four water-to-fuel ratios 


46/ J. A. Larrimore, footnote 14, pp. 185-210. 
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of the Cornell University Zero Power Reactor with sufficient 
accuracy to allow planning of further temperature dependence 
studies. 

The isothermal temperature coefficient of reactivity 
was measured for four water-to-fuel ratios with an accuracy 
which is comparable to that previously reported in the 
literature. The predicted and verified linear dependence 
ei 2 on FH,0 had not been reported previously. 

Since the experimental results are based on period 
measurements on a clean core, and since the calculations 
were done with a simple two group model, the disagreement 
between calculated and measured values of A and B indicates 


the need for a refinement in the calculations. 


6.2 TEMPERATURE DEPENDENCE OF THERMAL DISADVANTAGE FACTORS 
6.2.1 Summary 
The neutron density disadvantage factor OF and the 

experimentally measured disadvantage factor © exp are 
defined by the ratios ore = Wywt and oon = enyet where N 
and C are respectively the spatially averaged neutron 
density and foil count rate per unit volume; superscripts 
m and £f indicate respectively the moderator and fuel 
regions of the reactor. 

47 / 


MacVean's calculation of thermal lattice parameters— 


gives the temperature dependence of the disadvantage factor. 


47/ C. R. MacVean, footnote 2, p. 143. 
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Since the calculated changes are small, highly accurate 
experiments were required for comparison. Experimental 
disadvantage factors were measured for three water-to- 
fuel ratios over the temperature range 5-65°C. 

An integral foil counting method was used employing 
a 0.600+0.00l-inch diameter foil in the fuel and an 
appropriately shaped moderator foil. Five-mil thick 
foils of dysprosium (5.24% dysprosium in aluminum), 
manganese (P-metal), and vanadium were chosen because 
their cross sections are nearly 1/v. Two counting 
methods were used. One was to dissolve the manganese foils 
in nitric acid and gamma-count aliquots of the solutions 
in a well-type scintillation counter. The other method 
was to beta-count solid foils in a proportional flow 
counter. Results of the two methods agree within 
statistical errors. 

Experimental disadvantage factors were least-squares 
fitted to linear functions of temperature. Agreement 
between theory and experiment for the temperature depend- 
ence of the disadvantage factors is within the experimental 
error for all cores and foils, but the disadvantage factors 
themselves are a few per cent lower than the calculated 
disadvantage factors -- a result also found in the Brook- 


haven and Westinghouse work reported by Hence) 


48/ H. C. Honeck, footnote 10, p. 49. 
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622.25) Cone Misions 

Experimental determination of the thermal disadvantage 
factor was accomplished quite accurately using the technique 
of dissolving the foils and gamma-counting aliquots of the 
solutions. Beta-counting of the solid foils gave satis- 
factory results providing that measured values were corrected 
for the different geometrical counting efficiencies of 
different shapes of foils. The error for both methods was 
less than one per cent, an improvement on previously 
reported measurements. 

Theory and experiment disagree by a few per cent on 
the value of the disadvantage factor. This has been a 
common result, as mentioned above. Several factors may con- 
tribute to this disagreement: (1) flux perturbations due 
to finite thickness foils, (2) flux perturbations due to 
the lucite moderator foil support plate, (3) flux pertur- 
bations due to the use of activation rods with a shorter 
fuel region than normal fuel rods, (4) inexactly shaped 
moderator foils, and (5) a small error in the calculation 
of the disadvantage factor. The first two possibilities 
could be studied by using different thicknesses of foils 
and of moderator foil support plates. The third could be 
eliminated by using full length decappable rods, or 
minimized by placing foils only at the center of the fueled 


section of the activation rod. A more expensive die could 
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be made for punching the moderator foils to eliminate the 
fourth error source. Further study could be made on the 
calculational techniques for the disadvantage factor; 
however the estimated accuracy of 2% for the calculated 
values of the neutron density disadvantage factors by 


9 / 


MacVeant= indicate that the calculations probably would 
not change much with slight refinements to the code. 
Agreement between theory and experiment for the 
temperature dependence of the disadvantage factor is 
within the experimental error; however the experimental 
error in these results varies from about 20% to about 


90% because small differences were taken between relatively 


large numbers to obtain the temperature dependence. 


6.3 PROPOSALS FOR FURTHER STUDY 

When the perturbation theory calculations of the 
temperature coefficient are brought into agreement with 
the experiments, these same calculational techniques 
can be used to predict the effect of adding soluble 
poisons to the moderator. It would then be of interest 
to measure the temperature coefficient of reactivity 
as a function of the distributed poisons. 

As indicated earlier, it would be of great interest 


to investigate the effect that changing the foil thick- 


49 / Gs Re Macyveans, footnote 2, p. 111. 
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nesses would have on the disadvantage factor measurements. 
It might be of further interest to investigate the dis- 
advantage factor, and possibly its temperature variation, 
using resonance foils to study different energy regions. 
In conjunction with the above proposals, it would 
also be interesting to vary the thickness of the lucite 
moderator foil support plate and determine its effect on 
the measured value of the disadvantage factor. Another 
study could be made using polyethylene support plates. 
Further into the future, it might be of interest to 
measure these same parameters on cores which have different 
geometrical shapes. For example a square or long, thin 
lattice could be set up instead of a cylinder. The 
ability to properly predict the temperature coefficient 
for these geometries would be a measure of the accuracy 
of calculating the reflector contribution to the isothermal 


temperature coefficient of reactivity. 
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LATTICE PARAMETERS FOR THE CORNELL UNIVERSITY ZERO POWER 


REACTOR 


TAL ae 


CORNELL UNIVERSITY ZERO POWER REACTOR FUEL ELEMENT 


PARAMETERS 2/ 


Fuel Material 

Fuel Rod Diameter, in: 

Fuel Density, gm/cc: 

Fuel Enrichment, Weight Percent U 
Fuel Enrichment, Atom Percent yr??: 
Oxygen to Uranium Number Ratio: 
Fuel Molecular Weight; 

Cladding Material: 

Cladding Density, em/cc: 

Cladding Inner Diameter, in: 


Cladding Thickness, in: 


Fuel Element Diameter, in: 


a/ 


S. S. Berg, see footnote l. 
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235. 


Uranium Dioxide 
0.600+0.003 
L026 20.527 
220 0202.16 
22096090 
COL 
270.26(physical) 
6061-76 aluminum 
Zan 0 

0.610+0.002 
02028+0,003 
0.666+ .006 
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APPENDIX B 


ISOTHERMAL TEMPERATURE COEFFICIENT DATA REDUCTION 


Tables B.1 - B.4 give the following information in 
tabular form: date of the experiment, reference run 
number, number of fuel rods N, temperature 9 in °C, 
period T in seconds, excess reactivity oo aeeeS reference 
numbers of the two runs used to calculate Ap/A@, Ap, AQ, 
Ap /A® = a, the average temperature for the two runs 6, 
and the volumetric coefficient of thermal expansion of 
water at the average temperature “H,0° Figure B.l is 


a curve of PH o VS 98 obtained by differencing values of 
Z 


relative volume given in the Handbook of Physics and 


Chemistry 20/ 


oe ee ceo fe 


30/ Ge Us HOGeian, sl ootnote 27, p. 2143. 
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